Giant Radio Sources in View of the Dynamical Evolution of 
FRII-type Population. I. The Observational Data, and Basic 
Physical Parameters of Sources Derived from the Analytical 

Model 

by 

J. Machalski, K.T. Chyzy & M. Jamrozy 

Astronomical Observatory, Jagellonian University, ul. Orla 171, 
30-244 Cracow, Poland 
e-mail: machalsk@oa.uj.edu.pl jamrozy@oa.uj.edu.pl 

Abstract 

The time evolution of giant lobe-dominated radio galaxies (with 
projected linear size D > 1 Mpc if -ffo=50 kms _1 Mpc _1 and qo=0.5) 
is analysed on the basis of dynamical evolution of the entire FRII-type 
population. Two basic physical parameters, namely the jet power Qq 
and central density of the galaxy nucleus po are derived for a sample of 
giants with synchrotron ages reliably determined, and compared with 
the relevant parameters in a comparison sample of normal-size sources 
consisting of 3C, B2, and other sources. Having the apparent radio 
luminosity P and linear size D of each sample source, Qq and po are 
obtained by fitting the dynamical model of Kaiser et al. (1997). We 
find that: (i) there is not a unique factor governing the source size; 
they are old sources with temperate jet power (Qq) evolved in a rel- 
atively low-density environment (po). The size is dependent, in order 
of decreasing partial correlation coefficients, on age; then on Qq; next 
on po- (h) A self-similar expansion of the sources' cocoon seems to be 
feasible if the power supplied by the jets is a few orders of magnitude 
above the minimum-energy value. In other cases the expansion can 
only initially be self-similar; a departure from self-similarity for large 
and old sources is justified by observational data of giant sources, (iii) 
An apparent increase of the lowest internal pressure value observed 
within the largest sources' cocoon with redshift is obscured by the in- 
trinsic dependence of their size on age and the age on redshift, which 
hinders us from making definite conclusions about a cosmological evo- 
lution of intergalactic medium (IGM) pressure. 

Key words: galaxies: active - galaxies: evolution - galaxies: 
kinematics and dynamics 



1 Introduction 

Extragalactic radio sources, powered by twin jets resulting from nuclear en- 
ergy processes in the Active Galactic Nucleus (AGN), exhibit a very large 
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range of their linear size. The sizes of these powerful sources range from 
less than 10 2 pc (GPS: Giga Hertz-peaked spectrum), to 10 2 — 10 4 pc (CSS: 
compact steep spectrum), 10 4 — 10 6 pc (normal-size sources), up to greater 
than 10 6 pc = 1 Mpc ('giant' radio sources). One of the key problems of the 
evolution of extragalactic sources is whether and how different size sources 
are related. Is there a single evolutionary scheme governing the size evolution 
of radio sources, or do small and large sources evolve in a different way? 

From many years 'giant' 1 radio sources have been of special interest for 
several reasons. Their very large angular sizes give excellent opportunity 
for the study of source physics. They are also very useful to study the 
density and evolution of the intergalactic and intracluster environment (cf. 
Subrahmanyan and Saripalli 1993; Mack et al. 1998), as well as to verify 
the unification scheme for powerful radio sources (Barthel 1989; Urry and 
Padovani 1995). Finally, they can be used to constrain dynamical models of 
the source lifetime evolution (e.g. Kaiser and Alexander 1999). The general 
questions are: do the largest radio sources reach their extremal giant sizes 
due to (i) exceptional physical conditions in the intergalactic medium, (ii) 
extraordinary intrinsic properties of the AGN, or simply (iii) because they 
are extremely old? 

To answer these questions, in a number of papers attempts were made 
to recognize properties other than size which may differentiate giants from 
normal-size sources. The giant-source morphologies, energy density in the 
lobes and interaction with the intergalactic medium were studied by Sub- 
rahmanyan et al. (1996), who suggested that giant radio galaxies may be 
located in the lowest density regions, and 'may attained their large size as 
a result of restarting of their central engines in multiple phases of activity 
along roughly similar directions'. Also Mack et al. (1998), after a study of 
5 nearby giant radio galaxies, argued that those sources are so huge because 
of their low-density environment and not because of their old ages. A sim- 
ilar conclusion was drawn by Cotter (1998) after his study of a sample of 
high-redshift 7C giants. He found those sources to be not old and of similar 
kinetic powers of the jet as normal-size 3C sources. But the densities of their 
surrounding medium was found to be much lower than those around most 
of 3C sources (Rawlings and Saunders 1991). Ishwara-Chandra and Saikia 
(1999) compiled a sample of more than 50 known giant sources (many of 
them being of FRI-type morphology) and compared some of their properties 
with those of a complete sample of 3CR sources 'to investigate the evolution 
of giant sources, and test their consistency with the unified scheme for radio 
galaxies and quasars'. They concluded that the location of giants on the 

1 hereafter we use giant or giants instead of giant radio source(s) 
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power-linear size (P — D) diagram may suggest that the largest sources have 
evolved from the smaller sources. Finally, in the recent extensive study of 26 
giant galaxies by Schoenmakers et al. (2000), the authors argued that those 
galaxies 'are both old sources, in term of their spectral age, and are situated 
in a relatively low-density environment, but also that neither of these two 
properties are extreme. Therefore, their large size probably results from a 
combination of these properties'. 

From the above results, it is clear that the phenomenon of the giant radio 
sources is still open to further research. Therefore, in this paper we analyse 
whether observed properties of giant sources can be explained by a model 
of the dynamical evolution of classical double radio sources in cosmic time, 
and what factor (if there is a one) is primarily responsible for the giant size. 
Obviously such an analysis, based on the giant radio sources only, i.e. the 
sources with a strongly limited range of linear sizes, would not be reliable and 
even possible. Therefore, to solve the above problems we analyse the sample 
of radio sources comprising both the 'giant'-size and 'normal'-size sources, 
most of the latter with sizes from ~100 kpc to 1 Mpc. 

Most published analytical models of the dynamical evolution of extended 
powerful radio sources are based on the hydrodynamical self-similar expan- 
sion of the sources' cocoon caused by the interaction of light and supersonic 
jets with the ambient medium. Carvalho and O'Dea (2002) classify those 
models into three mutually exclusive types encoded as I, II, and III, and 
give an excellent description of their properties. However, we realize that 
only few of them deal with the source's energetics and the radio luminos- 
ity evolution with time which is crucial for e.g. the analysis of evolutionary 
tracks of sources on the luminosity-size (P-D) plane, i.e. the most sensitive 
characteristics of the dynamical models. 

Two type III models, published by Kaiser, Dennett-Thorpe and Alexan- 
der (1997) [hereafter referred to as KDA] and Blundell, Rawlings and Willott 
(1999) [hereafter BRW], are more sophisticated than models of type I and II. 
We apply the KDA model for its simplicity in comparison with the BRW one, 
and in spite of some objections about an application of the self-similar mod- 
els to large and old sources in which an internal pressure cannot be always 
above the external medium pressure (e.g. Hardcastle and Worrall 2000). 

The present analysis is confined to classical double radio sources with 
FRII-type (Fanaroff and Riley 1974) morphology. Synchrotron ages in a 
sample comprising both giant and normal-size sources are used to verify the 
dynamical time evolution of such sources predicted by the above analytical 
model. Basic physical parameters, i.e. the jet power Q , the central density of 
the galaxy nucleus po, the energy density and pressure in the lobes/cocoon 
(u c and p c ), and the total energy of the source E to t are derived with the 
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KDA model for each member of the sample to fit its estimated age, redshift, 
monochromatic radio luminosity, projected size, and axial ratio. Besides, 
the fitted values of some of these parameters are directly compared with 
the relevant values calculated straight from the data, e.g. the equipartition 
magnetic field strength, B eq , equipartition energy density in the source, u eq , 
and its total energy U eq . These 'observational' values and their errors are 
homogeneously calculated using the method outlined by Miley (1980). Next, 
the physical parameters derived for real radio sources are used to specify the 
conditions or circumstances under which the FRII-type radio sources can 
reach the largest observed linear sizes, to verify the applicability of the self- 
similar analytical models for the largest sources, and to search for an evidence 
after a cosmological evolution of the ambient pressure in the intergalactic 
medium (IGM). 

The observational data used and the physical parameters of the sample 
sources found directly from the data are given in Section 2. The application 
of the dynamical models is described in Section 3, while in Section 4 results of 
the modelling are presented in the form of statistically significant correlations 
between physical parameters of the sample sources derived from the fitting 
procedure and their apparent observational parameters. The results obtained 
are discussed in Section 5 and the conclusions are given in Section 6. 

2 Observational Data 
2.1 Selection Criteria 

Similarly to the approach of Ishwara-Chandra and Saikia, we have compiled 
a subsample of 18 giant sources and a comparison subsample of 54 normal- 
size sources. The selection criteria were as follows: (i) the sources have the 
FRII-type morphology, (ii) The existing radio maps enable a suitable deter- 
mination of their lateral extend, i.e. transversal to the source's axis, and (iii) 
their spectral age or the expansion velocity, determined with the same model 
of the energy losses, are available from the literature. The spectral ageing 
data calculated with the JP model (Jaffe and Perola 1973; cf. Sect. 2.4) for 
the giants are taken from the papers of Saripalli et al. (1994), Parma et 
al. (1996), Mack et al. (1998), Schoenmakers et al. (1998, 2000), Ishwara- 
Chandra and Saikia (1999), Lara et al. (2000), and Machalski and Jamrozy 
(2000). For the aim of this paper, i.e. for an observational verification of the 
dynamical time evolution of classical radio sources, especially the growth of 
their linear size with time predicted by the analytical models, the compari- 
son subsample has been chosen to comprise high-luminosity (high- and low- 
redshift), as well as low-luminosity normal-size sources. The high-redshift 
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(with z >0.5) and low-redshift (z <0.5) sets consist of 3C sources taken from 
the papers of Alexander and Leahy (1987), Leahy et al. (1989), Liu et al. 
(1992), and Guerra et al. (2000). All of them have P 178 > 10 25 WHz^sr" 1 
(other selection criteria are summarized in Liu et al.). The low-luminosity set 
comprises FRII-type sources with P 1A < 10 24,4 WHz~ 1 sr~ 1 (corresponding 
to P 178 < 10 25 WHz _1 sr _1 assuming a mean spectral index of 0.7 between 
178 and 1400 MHz). A limited number of such sources with spectral ages 
determined have been available from the papers of Klein et al. (1995) and 
Parma et al. (1999). 

2.2 Observational Parameters 

For all individual sample sources the following observational parameters are 
determined: the redshift z, the 1.4 GHz luminosity in WHz _1 sr _1 , the pro- 
jected linear size D in kpc, the cocoon's axial ratio AR = D/b and volume 
V Q in kpc 3 . The volume of the source, V Q , is calculated assuming a cylindrical 
geometry with the length D, and the base diameter b taken as the average 
of the full deconvolved widths of the two lobes. The latter are measured 
between 3a contours on a radio contour map half-way between the core and 
the hot spots or distinct extremities of the source. All these data for the 
giant-size and normal-size sources in our sample are given in columns 3-7 
of Table 1. The columns 8-9 of Table 1 give the reference papers to the 
radio map used to determine both D and AR for each sample source and its 
spectral age, respectively. 

2.3 Physical Parameters Derived Directly from the Data 

As mentioned in Section 1, the equipartition magnetic field B eq , the energy 
density u cq , and the total emitted energy U cq with their errors are calculated 
with the method outlined by Miley (1980). However, the Miley's assumption 
of pure power-law radio spectrum has been abandoned; the cocoon's radio 
spectrum has been determined by the least-square method fit of the simple 
analytic functions: y = a + bx + cexp(±x) or y = a + bx + ex 2 (where 
x = log z/[GHz], y = log S(u)) to the available flux densities S(u) weighted by 
their given error. The total luminosity of the cocoon has been then integrated 
between 10 MHz and 100 GHz using the above fitted spectrum with Hq = 
50 kms~ 1 Mpc~ 1 and g =0.5 2 . The values of -u cq , B cq , and t/ C q=w cq V r with 

2 An application of the most recent cosmological constants will change numerical values 
of dimensions, power, ambient density, etc., but not relations between observational and 
physical parameters of the sources. The applied constants provide an easier comparison of 
the derived physical parameters of sources with those found in a large number of previously 
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their estimated error, calculated for each od the sample sources with the 
assumption of a filling factor of unity and equipartition of energy between 
electrons and protons, are given in columns 3-5 of Table 2, respectively. 
Table 2 contains the adopted age and physical parameters of the sample 
sources derived either from the observational data (columns 3, 4, and 5) and 
from the analytical KDA model (columns 6, 7, 8, and 9). 

2.4 Spectral Age 

Determination of the age of the radio sources is crucial to constrain any 
dynamical model of their time evolution. An apparent age of sources can be 
estimated from the ratio of their total emitted energy, U cq , determined under 
the 'minimum energy' condition and the observed power 



and S(u) is the observed flux density at different frequencies. The resul- 
tant age of source, being rather its upper limit, is usually greater than the 
synchrotron age of relativistic particles commonly determined from the stan- 
dard spectral-ageing analysis (e.g. Alexander and Leahy 1987). However, 
the time-dependence of various energy losses suffered by the particles causes 
different parts of the lobes or cocoon to have different ages. Besides, the radi- 
ation losses (and thus the synchrotron age) depend on the history of particle 
injection, the distribution of the pitch angle, etc. described by the different 
synchrotron models: Kardashev-Pacholczyk (KP); Jaffe and Perola (JP); or 
continuous injection (CI), cf. Carilli et al.(1991) for a detailed description. 
Therefore, these models of different energy losses give different spectral age 
estimates, and a comparison of ages of the sample sources ought to be made 
within the same synchrotron model. Moreover, the synchrotron ages (t syQ ) 
usually differ from the dynamical age (£dyn) estimated from the ram-pressure 
arguments (cf. Begelman and Cioffi 1989; Lara et al. 2000; Schoenmakers et 



An attempt to minimize the discrepancies between spectral and dynami- 
cal ages has been undertaken by Kaiser (2000). His 3-dimensional model of 
the synchrotron emissivity of the cocoon traces the individual evolution of 
parts of the cocoon and provides, according to the author, a more accurate 
estimate for the age of a source. Its application to the lobes of Cygnus A 
gave a very good fit to their observed surface brightness distribution. How- 
ever, since an application of the above model is confined to the sources for 
which their lobes are reasonably resolved in the direction perpendicular to 

published papers 
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the jet axis - we cannot use it for our statistical approach to the dynamical 
evolution of giant radio sources and have to rely on results of the standard 
ageing analysis. 

Basing on a commonly accepted assumption about the proportionality of 
the spectral and dynamical ages, hereafter we assume td yn = 2t syn (e.g. Lara 
et al. 2000). This age (marked by t[Myr]) is given in column 2 of Table 2. 



3 Application of the KDA Model 

3.1 Source Dynamics 

The overall dynamics of a FRII-type source (precisely: its cocoon) described 
in the KDA model is based on the earlier self-similar model of Kaiser and 
Alexander (1997) [hereafter referred to as KA]. It is assumed that the radio 
structure is formed by two jets emanating from the AGN into a surrounding 
medium in two opposite directions, then terminating in strong shocks, and 
finally inflating the cocoon. A density distribution of the unperturbed ex- 
ternal gas is approximated by a power-law relation pa = Po(^/ao) , where 
d is the radial distance from the core of a source, po is the density at the 
core radius ao, and (3 is the exponent in this distribution [the simplified King 
(1972) model]. 

Half of the cocoon is approximated by a cylinder of length L } = D s /2 and 
axial ratio = AR/2, where D s is its total unprojected linear size. The 
cocoon expands along the jet axis driven by the hot spot pressure ph and 
in the perpendicular direction by the cocoon pressure p c . In the model the 
rate at which energy is transported along each jet (Qo) is constant during 
the source lifetime. The model predicts self-similar expansion of the cocoon 
and gives analytical formulae for the time evolution of various geometrical 
and physical parameters, e.g. the length of the jet [cf. equations (4) and (5) 
in KA]: 

Lj = DJ2 = ci ( — ^9 J t 3/(5 ^; (1) 



13 

Poa , 



and the cocoon pressure (cf. equation 34 in KA): 

where c\ is a dimensionless constant [equation (25) in KA], and T x is the 
adiabatic index of the unshocked medium surrounding the cocoon, and 7\c = 
Ph/pc is the pressure ratio. 
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However, the pressure ratio V^c = 4i?^, implied in the original KDA pa- 
per, has later been found to seriously overestimate the value of 7\c obtained 
in hydrodynamical simulations by Kaiser and Alexander (1999). Therefore, 
in our modelling procedure we use the empirical formula taken from Kaiser 
(2000): 

V hc = (2.14 - 0.52/3)^ 04 - a25/3 . (3) 



3.2 Source Energetics and Radio Power 

In our application of the KDA model we neglect thermal particles, hence the 
overall source dynamics is governed by the pressure in the cocoon in the form 
p c = (r c — l)(u c + Mb), where T c is the adiabatic index of the cocoon, u e and 
■ub are the energy densities of relativistic particles and the magnetic field, 
respectively. Both energy densities are a function of the source lifetime t. In 
particular, 

u B (t) oc B 2 {t) = const r a , (4) 

where a = (r B /r c )(4 + (3)/ (5 - (3). 

Since the time evolution of the pressure p c is known from the self-similar 
solution, then one can calculate the energy density in the cocoon at any 
specific age t: 

u c (t) = u e (t) + u B (t) = Pc (t)/(r c - 1), (5) 

and the total source energy: E tot (t) = u c (t)V c (t), where V c is the cocoon 
volume attained at the age t: 

V c (t)=2^\L i (t)f^fl^\ (6) 

Following KDA and Kaiser (2000) we can write: 

Etot = UcVc = 9[r c + (r c - 2 f)(nc/4)]-4-^ got - (7) 

Thus, the ratio of energy delivered by the twin jets and stored in the cocoon 
is: 

2Q t 9r c - 4 - p 9(r c -i) 

E tot 5-(3 + 4(5-/5) /hc ' 

i.e. for given T c and (3 values, this ratio is a function of the pressure ratio 
Phc only. For T c = 5/3 and (3 = 3/2 we have 
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2Q */£tat = 2.7 + 0.437\ c . (8) 

The radio power of the cocoon P u is calculated in the KDA model by 
splitting up the source into small volume elements and allowing them to 
evolve separately. The effects of adiabatic expansion, synchrotron losses, and 
inverse Compton scattering on the cosmic microwave background radiation 
are traced for these volume elements independently. The total radio emission 
at a fixed frequency v is then obtained by summing up the contribution 
from all such elements, resulting in an integral over time [equation (16) in 
KDA]. It depends on the source's age t and redshift z, the jet power Qq, the 
cocoon axial ratio R T , the exponent in the expected power-law distribution 
of relativistic particles p, and on the ratio r of the magnetic field energy 
to the energy of relativistic electrons and non-radiating particles (given in 
Sect. 3.3). The integral is not analytically solvable and, following the KDA 
we calculate it numerically. 

3.3 Fitting Procedure 

On the basis of the above model we aim to predict the specific physical 
parameters for all giants and normal-size sources in the sample at their es- 
timated (dynamical) age, i.e. Qo, po, u c , Pc, and E tot . This differs from the 
KDA approach, who on the base of available observational data, evaluated 
some general trends and made crude estimates of possible ranges of values 
attained by the model parameters. 

In order to derive the above parameters for our sources, all other free 
parameters of the model (r, p, T x , T c , Tb, ao, P), and the inclination angle 
of the jet axis to the observer's line-of-sight 9 have to be approximated. 
Following the KDA, we adopt their 'Case 3' where both the cocoon and 
magnetic field are 'cold' (r c = Tb = 5/3) and the adiabatic index of the 
jet material and external gas is also 5/3. For the initial ratio of the energy 
densities of the magnetic field and the particles we use r = ub/u c = (l+p)/4, 
with the exponent of the energy distribution p = 2.14. 

The core radius a is one of the most difficult model parameter to be 
set up. Even careful 2-D modelling of a distribution of radio emission for 
well known sources with quite regular structures can lead to values of ao 
discrepant with those predicted by X-ray observations, the only presently 
available method to determine the source environment (cf. an extensive 
discussion of this problem in Kaiser 2000). In our statistical approach we 
assume ao = 10 kpc for all sources, a conservative value between 2 kpc used 
by KDA and 50 kpc found by Wellman et al. (1997). In Section 5.1 we 
discuss the consequences of other possible values of this parameter. 
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We also use a constant value of (3 for all sample sources taking (3 = 
1.5 for further calculations. This is compatible with other estimates of this 
parameter (e.g. Daly 1995) although much flatter than (3 — 1.9 adopted 
in the original KDA paper on the basis of Canizares et al.'s (1987) paper 
who found that value to be typical for a galaxy at about 100 kpc from its 
centre. A flatter density profile should be more adequate for distances of a 
few hundreds of kpc. 

Another free parameter of the model, the orientation of the jet axis with 
respect to the observer's line-of-sight 9 = 90° is assumed for all giants and 
9 = 70° for other sources. This latter value is justified by the dominance 
of FRII-type radio galaxies in our sample. In view of the unified scheme 
for extragalactic radio sources, an average orientation angle (0rg) — 69° for 
radio galaxies only was determined by Barthel (1989). The apparent linear 
size D of a radio source then yields the model cocoon size: 

L > = 2lL < 9 » 

Having fixed all these free parameters of the model, we find the jet power 
Qo and the initial density of external medium p for each individual sample 
source by iterative solution of the system of two equation: (i) equation (1) 
equated to equation (9) for the jet length L s , and (ii) the integral for the 
luminosity of the cocoon P v [equation (16) in KDA; cf. Section 3.1.2] 
requiring the match of the solution to the observed values of D and Pi. 4 , 
respectively. The above fitting procedure proved to give always stable and 
unique solutions. Then, from equations (2) and (5) we calculate other model 
parameters: the cocoon's pressure p c and its energy density u c , and from 
equation (7) the total cocoon energy E tot . The resultant values of Qo, po and 
p c for the sample sources are given in columns 6, 7, and 8 of Table 2. 

In Table 3 we summarize meanings of the observational and model param- 
eters characterizing the source's cocoon and used in the present analysis. All 
dimensions are given in the SI units except the cocoons' length and volume 
which are given in kpc and kpc 3 , respectively. Two quantities present in the 
text do not appear in Table 3: the source (cocoon) unprojected linear size D s 
and unprojected volume V c which imply from the apparent (observed) size 
D, axial ratio AR, and assumed inclination of the sources' axis 9. 
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4 Results of the Modelling 



4.1 Jet Power Qq and Core Density po 

In the KDA model the jet power and ambient density are independent in 
accordance with a physical intuition. A distribution of these parameters 
derived for different sets of the sample sources on the log(<5o)~log(po) plane 
is shown in Figure la. One can realise in Figure la that giants are not fully 
separated from other sources, where (i) among the sources with a comparable 
jet power Q , giant sources have an average central density po smaller than a 
corresponding central density of normal-size sources, (ii) giants have at least 
ten times more powerful jets than much smaller low-luminosity sources of a 
comparable p . 

Moreover, for a number of sources in the sample the derived values of 
their fundamental parameters Qo and po are very close, while their ages are 
significantly different. Thus in view of the model assumptions, they may be 
considered as 'the same' source observed at different epochs of its lifetime. 
Such bunches of three to five sources (hereafter called 'clans') are indicated 
in Figure la with the large circles. These clans have appeared crucial in 
a comparison of the observing data and the model predictions, and in the 
analysis of the giant-source phenomenon. More detailed analysis of these 
clans and their evolution will be given in Paper II of this series. 

The KDA and BWR models predict that the luminosities of mature radio 
sources decrease with their age. Therefore, more distant sources fall below the 
flux-density limit of a sample sooner than nearer sources, and in any sample 
the high-redshift sources will be younger and more luminous than the low- 
redshift ones. A significant anticorrelation between Q Q and age t, expected 
as a consequence of the above effect (called 'youth-redshift degeneracy' in 
BWR), is shown in Figure lb. 

4.2 Cocoon's Energy Density u c and Total Emitted En- 
ergy E tot 

In the KDA model, the energetics of the radio source is governed by the 
jet power Q , adiabatic index r c , and the pressure ratio V\ IC . Since T c is 
assumed constant for all the sources and Q is constant for a given source, 
the energy of the cocoon u c is determined by the pressure p c attained by 
the cocoon at age t. As it decreases with time [cf. Equations (5), (3) and 
(2)] and the volume increases with time [Equation (6)], their product, i.e. 
the model total emitted energy, E to t, is the increasing function of time, and 
its value is a fraction of energy delivered by the jet since the source's birth, 
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Figure 1: Plots of the jet power Qo a) against central density of the core 
Po', b) against source age. The giants are indicated by crosses, high-redshift 
sources - diamonds, low-redshift sources - open circles, and low-luminosity 
sources - stars. The dotted lines mark a constant linear size predicted from 
equation (10). The numbers above some symbols indicate observed size of 
the marked source. 
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Figure 2: Plot of the energy density u c against total emitted energy E tot . 
The symbols indicating sources are the same as in Figure 1. Dotted lines 
mark the time axes for different values of Q . 



Qot [cf. Equation (7)]. A distribution of u c and .Etot parameters on the 
log(w c )-log(£'tot) plane is shown in Figure 2. The time axes, calculated with 
Equations (2), (5) and (7) for a constant jet power, are indicated by the 
dotted lines. 

In order to investigate the time evolution of giant radio sources on the 
basis of dynamical evolution of the entire FRII-type population, in the next 
Subsection we examine several correlations between the basic physical pa- 
rameters of the sample sources derived from the data and the above two 
models. In the result we realize that all statistical tendencies are similar in 
both models. Therefore, below we present these correlations between param- 
eters derived with the preferred KDA model only. 

4.3 Correlations between Observed and Model Param- 
eters 

Below we analyse the relations between principal observational and model 
parameters of the sample sources important for their time evolution. Most of 
these parameters are interdependent (for example: the linear size of a source 
is likely dependent both on its age and the ambient medium density), hence 
each parameter of sources in our sample correlates somehow with other pa- 
rameters. Therefore, in order to determine which correlation is the strongest, 
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we calculate the Pearson partial correlation coefficients between selected pa- 
rameters. For the reason that most correlations between different parameters 
seem to be a power law, all correlations are calculated between logarithms of 
the given parameters (for the sake of simplicity, the 'log' signs are omitted 
in all Tables showing the partial correlations). Hereafter r X y denotes the 
correlation coefficient between parameters X and Y , r X Y/w is the partial 
correlation coefficient between these parameters in the presence of a third 
one, (W), which can correlate with both X and Y, and Pxy/w is the prob- 
ability that the test pair X and Y is uncorrelated when W is held constant. 
Similarly, r X Y/vw, t xy/uvw, Pxy/vw, and Pxy/uvw are the correlation coef- 
ficients for the correlations involving four or five parameters, and the related 
probabilities, respectively. 

A strong correlation between linear size and spectral age of 3C radio 
sources was already noted by Alexander and Leahy (1987) and confirmed by 
Liu et al. (1992). This correlation in our samples is shown in Figure 3. The 
giant sources do not show any tendency to a faster expansion velocity than 
that for the normal-size sources. The same conclusion has been made by 
Schoenmakers et al. (2000). However, two other aspects are worth empha- 
sizing: (i) There are four high-redshift giants which are much younger than 
the low-redshift giants. Two of them are quasars. It seems that they might 
grow so large under some exceptional conditions, (ii) The D — t relation for 
the low-luminosity sources (mostly B2) follows the same slope of the correla- 
tion as that for other sources, but low-luminosity sources are definitely much 
smaller indicating a dependence of the size and expansion velocity on the 
source luminosity. 

The partial correlation coefficients between the size D and t, Qq, and po 
together with the related probabilities of their chance correlation are given 
in Table 4. 

In view of the dynamical model applied and as a result of the above statistical 
correlations, we see that the linear size of a source strongly depends on both 
its age and the jet power, where the correlation with age is the strongest. 
However, the size also anti-correlates with central density of the core. That 
anticorrelation seems to be a weaker than the correlations with Qq and t and 
become well pronounced only when all three remaining parameters (Qq, t 
and z) are kept constant. 

Fitting a surface to the values of D over the t-Qo plane, we find 

D(t,Q ) Kt im±om Q° 31±002 . (10) 

The above relation well illustrates the influence of the jet power on the source 
(cocoon) expansion velocity, i.e. its length at a given age. However, the 
source expansion velocity ought to depend also on the external environment 
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Figure 3: Plot of the linear size D against source age t. The symbols indicat- 
ing sources are the same as in Figure 1. The solid lines indicate the implied 
expansion velocities in the speed of light c units. The dashed line shows the 
model predicted D{t) oc Lj(t) relation resulting from Equation (1) for (3—1.5. 

conditions. Indeed, the significant partial correlation coefficient for the D—p 
correlation in Table 4 confirms this effect. 

A greater than one exponent of the age in Equation (10) may suggest a 
statistical acceleration of the expansion velocities with age which contradicts 
with the deceleration implied by Equation (1) for ,3=1.5 (shown with the 
dashed line in Figure 3). However, the observed luminosities of the sample 
sources demand different jet powers (cf. see below). A linear regression of 
their apparent D values transformed to a constant reference value of Qo with 
Equation (10) on the age axis gives D(t) oc t L05±0 - 05 . This is somehow a 
surprised result and there would be several explanations of it: (i) (3 > 2, 
(ii) non-constant Qo during the lifetime of sources, (iii) not representative 
sample of sources, or (iv) real possibility of an acceleration of the expansion 
speed, at least for sources evolving in a specific environmental conditions. 
We will return to this possibility in Paper II. 

The partial correlation coefficients calculated for the correlation between 
the luminosity Pi. 4 and t, Qo and 1+z, as well as the related probabilities of 
their chance correlations are given in Table 5. 

Table 5 shows that the strongest correlation is between the source apparent 
luminosity and its jet power, although the anticorrelation between the lumi- 
nosity and age is also significant. The model values of Qo vs. the observed 
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Figure 4: Plot of the model jet power Q against the observed 1.4 GHz 
luminosity of the sample sources. 



1.4-GHz luminosities of the sample sources are plotted in Figure 4. 

The giant sources seemingly have either a higher jet powers or lower radio 
luminosities than the normal-size sample sources. According to the model 
applied, the source luminosity (at a constant Q ) decreases with time, thus 
Figure 4 confirms that the giants, statistically older than normal-size sources 
in the sample (cf. Table 2), are less luminous than younger comparison 
sources with a similar Q . 

Fitting a surface to the values of P 1A over the t-Q plane, we find 



PlA(t, Qo) OC t 



-0.61i0.08/n0.99i0. 05 



Qo 



An anticorrelation between apparent luminosity and age of matured sources 
has been predicted by the KDA and BWR analytical models. The data of the 
sample sources allow to verify those predictions. This anticorrelation in our 
sample is shown in Figure 5 where the 1.4-GHz luminosities are transformed 
to the constant jet power of 10 38 W according to the surface fit given above. 
There is an evident lack of powerful old sources concordant with predictions 
of the above two models. However, the luminosity of giant sources seems to 
decrease faster with respect to that predicted by the KDA model and very 
likely is connected with a departure from the self-similarity of the cocoon 
expansion in that model. Therefore our data on giant sources support rather 
the evolutionary predictions of the BWR model. 

The first time (to our knowledge) axial ratios of giant sources were anal- 
ysed and compared with those of smaller FRII-type sources by Subrah- 
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Figure 5: Plot of the 1.4 GHz luminosity P L4 transformed to a constant jet 
power of 10 38 W against the source age t. The symbols indicating sources 
are the same as in Figure 1. The dashed curve shows the P — t relation 
for Qo=10 38 W resulting from the iterative solution of the relevant KDA 
equations; cf. Section 3.3. 



manyan et al. (1996), who found no difference between the axial ratios of 
eight giants and eight 3C sources with a median size of about 400 kpc. The 
authors did not specify which 3C sources were considered, but since all giant 
and normal sources were of comparable powers and at comparable redshifts, 
we assume they might be of similar ages, so the dependence of AR on time 
could not be visible. 

In the BRW model the axial ratio of an individual source steadily increases 
throughout its lifetime. Moreover, that model implies a dependence of the 
AR on the jet power Q . The latter dependence was probably reflected by an 
apparent correlation between AR and the 178-MHz luminosity of 3C sources 
noted by Leahy and Williams (1984). Taking into account the unavoidable 
anticorrelation between Q and age in any sample of sources (cf. Section 4.1), 
in Table 6 we have calculated the partial correlation coefficients and the 
related probabilities of chance correlations between AR and t, AR and Qo, 
and AR and p when relevant combinations of the parameters t, Qo, po, and 
1+z are kept constant. 

Table 6 shows statistically significant correlations between the axial ratio 
and the source's age, as well as the axial ratio and the jet power. Fitting a 
surface to the values of AR over the Qo-t plane (where Q is in watts and t 
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Figure 6: Plot of the source axial ratio transformed to a constant jet power 
of 10 39 W against the age t. The dashed line indicates the linear regression 
on the age axis. 



in Myr), we found 

AR{t,Q ) <xt°- 23±om Q° 12±002 . (11) 

Indeed, our statistical data strongly support the implication of the BRW 
model of the dependence of AR on Q . A consequence of this effect for the 
expansion speed of the cocoon is pointed out in the next subsection. Using 
the above relation, we transform the apparent AR values from Table 1 to 
a reference jet power of 10 39 W. The relation between the transformed axial 
ratio and age of the sample sources with the regression line on the time axis 
is shown in Figure 6. 

This statistical correlation between the cocoon's axial ratio and age im- 
plies a time evolution of the ratio of the pressure in the head of the jet to 
the cocoon pressure. Indeed, substitution of Equation (11) into Equation (3) 
(for (3—1.5 with its assumed uncertainty of ±0.4) gives: 

n c (t,go)oct a38±0 - 06 gr i±a04 - (12) 

which violates the model assumption of self-similar expansion of the cocoon. 
The consequence of this will be analysed in Paper II. 

Our statistical analysis confirms the significant anticorrelation between 
the cocoon pressure and the size of sources expected from Equation (2), 
but reveals also a significant correlation between this pressure and redshift. 
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Figure 7: a) Cocoon pressure against size of the sample sources; b) the 
same pressure transformed to the reference size of 1 Mpc and age of 100 
Myr against redshift. The symbols indicating sources within different sets 
are the same as in Figure 1 and 2. The dashed line in a and b indicates the 
pressumed IGM presure evolution J9igm oc (1 + z) 5 . 
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Besides, the size strongly correlates with age and anticorrelates with redshift, 
so we have calculated the partial correlation between all these parameters. 
The Pearson partial correlation coefficients between p c , D, t, and 1+z are 
given in Table 7. 

Table 7 shows that: (i) The correlation coefficients for the direct (t X y) 
and partial (txy/v, t xy/uv) correlations between the cocoon pressure and 
its size, and the cocoon pressure and the source age are very similar. This is 
obvious in view of the very high size-age correlation shown in Table 4. (ii) 
Both direct correlations are weakened, though still remain very significant, if 
the age and redshift (in the p c — D correlation) and the size and redshift (in 
the p c — t correlation) are kept constant. Formally, the strongest partial cor- 
relation is found between the cocoon pressure and size (in fact: the cocoon's 
volume), (iii) the strong direct correlation between the cocoon pressure and 
redshift is seriously weakened if the source size and age are kept constant. 
This correlation, shown in Figure 7a, is very important for studies of physi- 
cal conditions in the intergalactic medium (IGM). The dashed line indicates 
the expected electron pressure in an adiabatically expanding Universe in the 
form Pigm = Pigm(1 + z ) 5 with Pi GM = 2 • 10~ 15 Nm" 2 (cf. Subrahmanyan 
and Saripalli 1993). 

Fitting a surface to the values of p c over the D — t plane we find 

p c (D,t)ocD-°- 30±0 - 19 t- 1 - 78±0 - 17 . 

Therefore, one can transform the cocoon's pressure values to a reference size 
and age. The plot of cocoon pressures transformed to D=l Mpc and t=100 
Myr versus 1+z is shown in Figure 7b. This Figure illustrate how p c of 
the high-redshift sources would decrease if they had evolved into the above 
size and age, and emphasize that the strong direct p c —{l+z) correlation in 
Figure 7a is, in fact, due to the much stronger correlations p c — D and p c — t 
(cf. Table 7). However, this is purely statistical result and the question about 
conditions under which distant radio sources would reach very large size and 
old age remains open. The aspect of a cosmological evolution of the IGM is 
discussed in Section 5.4. 



5 Discussion of the Results 

5.1 Influence of Fixed Model Parameters on the Model 
Predictions 

The basic physical parameters of the radio sources derived in this paper 
with the aim of the KDA model, i.e. jet power Q , central density po, and 
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cocoon pressure p c , are in principle dependent on the assumed central core 
radius, the exponent in the external gas distribution, the adiabatic indices 
of electrons and magnetic field, as well as on the orientation of the jet axis 
towards the observer. In application of the KDA model we have assumed the 
same values of these parameters for all sample sources. This assumption can 
only be valid in our statistical analysis of the evolutionary trends in the whole 
FRII-type population but not for individual sample sources. In particular, 
we have adopted do = 10 kpc and (3 = 1.5, respectively. Taking a lower core 
radius, e.g. a = 2 kpc, and keeping (5 — 1.5 will result in increase of the 
model density of the core po by roughly one order of magnitude, while other 
model parameters will not be changed. Conversely, a lower density gradient, 
e.g. (3 = 1.1, will lower p approximately 1.5 ~ 6 times. In this case however, 
the jet power will be increased by a few percent and accordingly the pressure 
and energy density in the cocoon will be changed. Solving the equations 
in Section 3 for a few sets of the model free parameters (a , (3, 0) we find 
that the presented and discussed correlations between the observational and 
model parameters are not changed; all the statistical trends are preserved 
although the values of the model parameters (especially p ) are changed 
quantitatively. 

A relativistic equation of state for the magnetic field does not change our 
results significantly, unless the KDA 'Case 1' (T c = T B = 4/3) is adopted. 
However, this case, i.e. when both the cocoon and the magnetic field energy 
have a relativistic equation of state, is unlikely for our sample sources. There- 
fore, we argue that the results discussed below are not significantly biased 
by selection of a particular set of the model parameters. 

5.2 Cause of Extremal Linear Size 

In view of the KDA and BRW analytical models of dynamical evolution of 
FRII-type radio sources, many such sources can evolve into a stage charac- 
terized by a linear size exceeding 1 Mpc. Access to this stage depends on 
a number of the model parameters: jet power Q , its Lorentz factor 7j et , 
the adiabatic indices of the cocoon material and magnetic field, T c and Tb, 
respectively, as well as the core radius ao, the external gas density po, and 
the exponent of its distribution (3. For a given set of these parameters, the 
model allows us to determine whether an evolving source will reach the size 
of 1 Mpc, and if so, at what age. From this point of view, giant sources 
should be the oldest ones. 

In this paper (Section 4.3) we have confronted the KDA model predictions 
with the observational data on giant-size and normal-size FRII sources. Our 
statistical analysis strongly suggests that there is not a single evolutionary 
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scheme governing the size development. An old age or a low external density 
alone is insufficient to assure extremely large linear extent of a source, both 
are necessary together with a suitable power driven from AGN by the highly 
relativistic jets. The Pearson partial correlation analysis indicates that the 
dependence of linear size on each of these three parameters is statistically 
significant. Ordering these correlations by decreasing partial correlation co- 
efficients, we find that the size is dependent on age, then on Q ; next on 
Po- 

About 83 % of giants in our sample possess a projected linear size over 1 
Mpc owing to statistically old age, low or moderate density of the external 
medium, and high enough power of their jets. The remaining 17 % (3 sources) 
are high luminosity sources at redshifts z > 0.4 ~ 0.5 and ages from 15 to 
25 Myr which are typical for normal-size sources. Two of them are quasars. 
The jet power of these sources is high enough to compensate for a higher 
ram pressure in a denser surrounding environment and higher energy losses 
during the cocoon expansion. The jet power of giants is not extreme, so 
several FRII-type sources having that Qo can potentially achieve very large 
size after a suitably long time. According to our results (cf. Figure la), 
these potential giants should have Qo > 10 375 W and be situated in an 
environment with po < 10~ 23 kgm~ 3 . 

The above scenario is not caused by a selection effect. We show that 
there are low-luminosity sources which lie in the parts of log Q — log p 
plane (Figure la) completely avoided by giant sources. They diverge from 
giants and even normal-size powerful sources by having low jet power Q < 
10 37 ' 5 W and total energy E tot < 10 52 - 5 J. Thus we suggest that they never 
reach the giant size. It is worth noting that some of them already have 
an age comparable to that of typical giants, in accordance with the model 
expectations. 

5.3 Jet Power and Energy Budget 

The ratio of total energy supplied by the twin jets during the lifetime of a 
source (2Q t) and its energy stored in the cocoon, derived from the data un- 
der the assumption of energy equipartition, U cq = u cq V , allow another test 
of the dynamical model predictions. The ratio of 2Q t/U cq for the sample 
sources (given in column 9 of Table 2) is plotted in Figure 8 vs. the cocoon 
axial ratio AR. The uncertainties in both values are marked by error bars 
on some data points. The solid curve indicates the model prediction from 
Equation (8), while the dashed curve shows the best fit to the data. The 
observed trend fully corresponds to the model prediction. However, the de- 
rived values of 2Q t/U eq , i.e. the reciprocal of the efficiency factor by which 
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Figure 8: Ratio of the total energy supplied by twin jets and energy stored 
in the cocoon against its axial ratio AR. The large uncertainties in both 
parameters are marked by error bars for a few sources only for clarity. Here 
the giant sources are marked with filled circles. The solid curve indicates the 
model prediction from Equation (8); the dashed curve shows the best fit to 
the weighted data. 

the kinematic energy of the jets is converted into radiation, is much higher 
than ~2, a value usually assumed in a number of papers. 

In a number of studies of giant radio sources (e.g. Parma et al. 1996; 
Schoenmakers et al. 1998) the authors followed the approach of Rawlings 
and Saunders (1991) and assumed a fraction of the total jet energy wasted 
for adiabatic expansion of the cocoon to be about 0.5 and used it to estimate 
the jet power Q for sources with known age (almost always from spectral 
ageing analysis). In the KDA model the energy stored in the source (cocoon) 
is 

E tot » J {Q dt - (p c d[V c (t)} + Ph d[V h (t)})} 

where p c dV c + ph dVh is the work done to expand the cocoon, and p^ and 
14 are the hotspot pressure and volume, respectively. If Vh is neglected, the 
expansion work will be 0.5 Qot] if not, it is dependent on the pressure ratio 
Phc as shown in Section 3.2. 

The values of 2Q t/U eq in our sample (Figure 8; Table 2) vary from about 
2 to more than 10 due to the high jet power required for the adiabatic ex- 
pansion of the cocoon. The Q values derived here from the KDA model can 
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be compared with the relevant values estimated by Wan, Daly and Guerra 
(2000) [WDG] for 22 3C sources included also in our sample (21 with z>0.5 
+ CygA). Recalculating their values for Hq=50 kms _1 Mpc _1 we found the 
WDG estimates approximately 2.5 times lower than the KDA values. Ex- 
actly the similar ratio (1.7-S-5.5 depending on the value of 7\c) characterize 
the jet powers of 9 3C sources, common with our sample sources, determined 
by Rawlings and Saunders (1991). The explanation of this ratio is straight- 
forward. The values of Qq estimated in those papers were based on the 
ram pressure considerations in the overpressured source model A of Scheuer 
(1974) and its further modifications (e.g. Begelman and Cioffi 1989; Loken 
et al. 1992; Nath 1995). They all are self-similar models of the Carvaldo 
and O'Dea type I which describe the source dynamics only. If the source 
energetics and especially the energy losses are properly taken into account 
(the type III models; e.g. KDA, BRW), the significantly higher values of Qo 
are implied. 

The data in our sample show a dependence of the energy ratio 2Q t/U cq on 
the cocoon axial ratio AR, and imply an increase of the fraction of jet energy 
spent on the adiabatic expansion of the cocoon volume in time. However, 
the data also suggest that for a constant AR (i.e. a given geometry of the 
cocoon), giants tend to have a smaller ratio of 2Q t/U eq than normal-size 
sources which means less energy of the jets converted into adiabatic expansion 
of the cocoon. This may indicate a lower pressure of the external medium 
surrounding the giant sources than that around smaller ones. 

5.4 External Pressure of the Surrounding Medium and 
its Evolution 

A non-relativistic uniform intergalactic medium (IGM) in thermal equilib- 
rium filling an adiabatically expanding Universe should have an electron pres- 
sure evolving with redshift Pigm = Pigm(1 + - 2 ) 5 - The advancing hotspots of 
FRII-type radio sources are probably confined by ram pressure of the IGM. 
Giant sources, with their lobes extended far outside typical galaxy halo, have 
the lowest values of p c and may be useful for determining the upper limit of 
Pigm- Using a small sample of giant sources, Subrahmanyan and Saripalli 
(1993) limited its local value to p^ GM w (0.5 + 2) • 10" 15 Nm" 2 . A further 
study was undertaken by Cotter (1998) who, using a larger sample of 7C 
giants with sources out to redshift of ~0.9, confirmed a strong dependence 
of the lowest p c on redshift in agreement with a (1 + z) 5 relation. This obser- 
vational result has been critically discussed by Schoenmakers et al. (2000), 
who have considered possible selection effects in Cotter's analysis (including 
the Malmquist bias), and concluded that there was not evidence in their 
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own sample for a cosmological evolution of Pigm- However, they also state 
that this hypothesis cannot be rejected until some low-pressure high-redshift 
sources are found. 

In all the above analyses the age of sources was not considered. The very 
significant correlation between p c and D or t, shown in Section 4.3, strongly 
suggests that the intrinsic dependences of size on age as well as of age on 
redshift and not the Malmquist bias is mainly responsible for the apparent 
correlation between p c and 1+z. Nevertheless, we agree with Schoenmakers 
et al.'s conclusion that until giant sources with internal pressures in their 
lobes p c < 2 ■ 1CT 15 N mr 2 at redshifts of at least 0.6 0.8 are not discovered, 
the IGM pressure evolution in the form Pigm oc (1 + z) 5 cannot be rejected. 

Most giants in our sample, except four high-redshift ones, reveal the low- 
est pressure in their cocoons. Sharing Subrahmanyan and Saripalli's argu- 
ments, we can expect those cocoon pressures are indicative of an upper limit 
to the present-day external pressure of the IGM, Pi GM . Taking into account 
the lowest values of p c (cf. Figure 7a), we found Pi GM < 2 • 10~ 15 Nm~ 2 
in accordance with their value. It is worth emphasizing that the above re- 
sults are obtained from the analytical model assuming energy equipartition 
in the initial ratio of the energy densities of the magnetic field and the par- 
ticles (cf. Section 3.1). This may not be the case in every part of the 
source (cocoon). Hardcastle and Worrall (2000) estimated gas pressures in 
the X-ray-emitting medium around normal-size 3CRR FRII radio galaxies 
and quasars, and found that, with few exceptions, the minimum pressures in 
their lobes, determined under equipartition conditions, were well below the 
environment pressures measured from old ROSAT observations. Therefore, 
they have argued that there must be an additional contribution to the inter- 
nal pressure in lobes of those sources likely including pressure from protons, 
magnetic fields exceeding their minimum-energy values, or non-uniform fill- 
ing factors. Nevertheless, the diffuse lobes of giants, extending farther from a 
host galaxy than the typical radius of high-density X-ray-emitting gas, may 
be in equilibrium with an ambient medium the emissivity of which is not 
directly detectable. 

6 Conclusions 

In this paper we confront the analytical KDA model predictions with the 
observational data of 'giant' and normal size FRII- type radio sources. From 
our analysis we can conclude as follows: 

(1) Giant sources do not form a separate class of radio sources, and do 
not reach their extremal sizes exclusively due to some exceptional physical 
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conditions of the external medium. The size is dependent, in order of de- 
creasing partial correlation coefficients, on age; then on the jet power Q ; 
next on the central core density po- 

(2) Giants possess the lowest equipartition magnetic field strength and 
energy density of their cocoons making their detection difficult in synchrotron 
emission. However, their accumulated total energy is the highest among all 
sources and exceeds 3 • 1 52 W. 

(3) Our data confirm the conclusion drawn by Blundell et al. (1999) that 
throughout the lifetime of an individual source its axial ratio can steadily 
increase, thus its expansion cannot be self-similar all the time. A self-similar 
expansion seems to be feasible if the power supplied by the jets is a few orders 
of magnitude above the minimum-energy value. In other cases the expansion 
can only initially be self-similar; a departure from self-similarity for large and 
old sources is justified by observations of giant sources. 

(4) The apparent increase of the lowest cocoon pressures (observed in the 
largest sources) with redshift is mainly caused by the intrinsic dependence of 
their age on redshift and dominates over a bias by possible selection effects. 
However, a cosmological evolution of the IGM cannot be rejected until giant 
sources with internal pressures in their lobes less than 2 • 10~ 15 Nm 2 at high 
redshifts are not discovered. 
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Observational data and physical parameters of the sample sources 



Table 1 gives the observational data of the sample sources. Their radio 
luminosity, size and volume are calculated with H o =50 kms _1 Mpc _1 and 
g =0.5. The entries are: 

Column 1 : IAU- format name of source 

Column 2: Other name 

Column 3: Redshift 

Column 4: Logarithm of 1.4 GHz luminosity 

Column 5: Projected linear size 

Column 6: Source axial ratio 

Column 7: Logarithm of source volume 

Column 8 and 9: References to the radio maps and spectral ageing anal- 
ysis, respectively 

Table 2 contains the adopted age and physical parameters of the sample 
sources derived either from the observational data (columns 3, 4, and 5) and 
from the analytical KDA model (columns 6, 7, 8, and 9). The entries are: 

Column 2: Age of source 

Column 3: Logarithm of the equipartition energy density 
Column 4: Equipartition magnetic field strength 
Column 5: Logarithm of the source total energy 
Column 6: Logarithm of the jet power 

Column 7: Logarithm of the initial density of the external medium at the 
core radius a 

Column 8: Logarithm of the cocoon pressure 

Column 9: Ratio between the total energy of twin jets and the source 
energy. 
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Table 1: Data of the sample sources 
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OA 
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5.47±0.07 


21 


20 


1723+510 


3C356 


1.079 


26 


.96 


643±13 


7.9±1.0 


6.55±0.10 


9 


9 


1726+318 


3C357 


0.1664 


25 


.43 


395±10 


3.0±0.6 


6.76±0.16 


21,3 


20 


1957+405 


CygA 


0.0564 


27 


.25 


185±3 


3.8±0.5 
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Table 2: Age and physical parameters 



Source 
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[Myr] 
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-13.63 + 0. 


27 





.16+0 


.05 


53 


.24+0 


.36 


38. 


.26 


-23 


.21 


-13 


,73 


8, 


.9+5.7 


2043+749 


64+11 


-13.57 + 0. 


10 





.17+0 


.02 


53 


.04+0 


.23 


38. 


14 


-24 


.06 


-13 


,74 


5. 


,1+1.9 


NORMAL 




































0154+286 


13±2 


-11.55 + 0. 


16 


1 


.74+0 


.32 


53 


.46+0. 


.32 


39 


,43 


-22 


.98 


-11. 


,72 


7, 


,7+5.0 


0229+341 


11±2 


-11.31 + 0. 


19 


2 


.31+0 


.43 


53 


.65+0 


.35 


39. 


,69 


-22 


.91 


-11. 


,33 


7, 


,8+5.4 


0231+313 


4.6+0.3 


-10.71 ± 0. 


16 


4 


.58+0 


.85 


53 


.62+0 


.29 


39 


.64 


-23 


.55 


-11, 


,03 


3, 


,2+1.9 


0404+428 


18+3 


-11.91 ± 0. 


17 


1 


.15+0 


.23 


53 


.08+0 


.30 


39. 


.02 


-22 


.96 


-12. 


,01 


10. 


.2+5.5 


0610+260 


45+4 


-11.72 + 0. 


18 


1 


.43+0. 


.30 


53 


.47+0 


.36 


39. 


.00 


-22 


.20 


-11. 


,54 


10. 


,2+7.6 


0640+233 


60+10 


-12.34 + 0. 


16 





.70+0 


.13 


53 


.03+0 


.32 


38. 


,33 


-22 


.66 


-12, 


,30 


7, 


.9+4.5 


0642+214 


14+5 


-11.71 + 0. 


22 


1 


.44+0 


.37 


52 


.58+0 


.34 


38. 


,87 


-22 


.71 


-11, 


,23 


17. 


,9+7.9 


0710+118 


35+5 


-11.54 + 0. 


16 


1 


.76+0 


.32 


53 


.57+0 


.34 


39 


.09 


—22 


.31 


-11, 


,52 


7, 


.8+5.0 


0806+426 


7.0+1.2 


-10.46 + 0. 


12 


6 


.13+0 


.82 


53 


.22+0 


.33 


39 


.27 


-22 


.56 


-10. 


,57 


5, 


,0+2.3 


0828+324 


59+9 


-13.27 + 0. 


18 





.24+0 


.05 


51 


.90+0 


.27 


37. 


.20 


-23 


.60 


-13 


,24 


7, 


.8+3.6 


0908+376 


28+5 


-12.33 + 0. 


16 





.71+0 


.13 


51 


.38+0 


.23 


36 


,83 


-23 


.36 


-12. 


,39 


5, 


,5+2.0 


0958+290 


22+4 


-12.11 + 0. 


15 





.92+0 


.16 


52 


.94+0 


.30 


38. 


,53 


-23 


.31 


-12. 


,32 


5, 


,7+2.9 


1008+467 


2.8+0.3 


-10.28 + 0. 


16 


7. 


.48+1 


.39 


53 


.20+0 


.35 


39 


,66 


-23 


.21 


-10. 


,36 


5. 


.5+4.0 


1012+488 


38+6 


-12.83 + 0. 


13 





.40+0 


.06 


53 


.40+0 


.23 


38. 


,81 


-23 


.93 


-12. 


,80 


6. 


,7+2.6 


1030+585 


14+2 


-11.56 + 0. 


16 


1 


.72+0 


.32 


53 


.01+0 


.33 


38. 


,98 


-22 


.84 


-11. 


,62 


8, 


,9+5.5 


1056+432 


3.2+0.3 


-10.72 + 0. 


14 


4 


.55+0 


.71 


52 


.75+0 


.31 


39 


,18 


-23 


.44 


-10. 


,81 


5. 


.7+3.5 


1100+772 


32+6 


-11.92 + 0. 


16 


1 


.14+0 


.21 


52 


.76+0 


.38 


38. 


,24 


-22 


.87 


-11, 


.95 


6. 


,5+4.6 


1111+408 


3.1+0.2 


-10.75 + 0. 


15 


4 


.40+0 


.76 


52 


.94+0 


.24 


39 


,20 


-23 


.79 


-10. 


,97 


3 


.8+1.9 


1113+295 


19+3 


-12.42 + 0. 


16 





.64+0. 


.12 


51 


.25+0 


.26 


36 


,84 


-23 


.84 


-12, 


,51 


5. 


,0+2.2 


1140+223 


1.7+0.6 


-9.87 + 0.11 


12.0+1 


.50 


52 


.80+0 


.24 


39 


,38 


-23 


.27 


-10. 


,02 


4, 


,4+1.0 


1141+354 


3.4+0.8 


-10.69 + 0.16 


4.68+0 


.87 


52 


.71+0 


.26 


39 


,03 


-23 


.44 


-10. 


,85 


4, 


,8+2.4 


1142+318 


22+5 


-11.66 + 0.16 


1.54+0 


.29 


53 


.69+0 


.37 


39 


,64 


-22 


.52 


-11, 


,57 


13 


,0+8.1 


1143+500 


1.2+0.5 


-9.67 + 0.16 


15.1+2 


.80 


52 


.22+0 


.23 


39 


,32 


-22 


.65 


-9. 


.58 


10 


,0+4.2 
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Source 


t 




Igv 


cq 








B eq 










IgQo 


IgPo 




k 


5Pc 




2Q t 




[Myr] 


[Jm 


-3 


] 






[nT] 






[J] 




[W] 


[kgm 3 ] 


[Nm 


- 2 ] 




C/cq 


1147+130 


12+4 


-11 


.10 


± 


0. 


.14 


2. 


.93+0. 


,48 


53 


.55+0. 


.26 


39 


.51 


-22. 


.60 


-11 


.16 


7 


.7+2.4 


1157+732 


13+3 


-11 


.17 


± 


0. 


16 


2. 


.70+0. 


.49 


53 


.77+0. 


.27 


39 


.73 


-22. 


.73 


-11 


.20 


7. 


.8+3.2 


1206+439 


3.0+0.3 


-10 


.34 


± 


0. 


.12 


7, 


.00+1. 


.00 


52 


.98+0. 


.20 


39 


.52 


-23 


.04 


-10 


.32 


7 


.1+2.6 


1216+507 


50+15 


-13 


.11 


± 


0. 


18 


0. 


.29+0 


,06 


52 


.75+0. 


.28 


37 


.98 


-23 


.73 


-13 


.29 


5 


.8+1.9 


1218+339 


3.9+1.1 


-10 


.35 


± 


0. 


.17 


6. 


.90+1. 


,35 


53 


.26+0 


.30 


39 


.60 


-23 


,00 


-10 


.41 


5 


.7+2.4 


1221+423 


18+8 


-11 


.94 


± 


0. 


.16 


1. 


.12+0. 


.21 


53 


.41+0. 


,36 


39 


.25 


-23 


,26 


-11 


.90 


8. 


.1+3.1 


1241+166 


4.2+0.4 


-11 


.40 


± 


0. 


.16 


2. 


,07+0. 


,38 


52 


.73+0. 


,30 


39 


.01 


-24. 


,21 


-11 


.44 


5 


.5+3.3 


1254+476 


2.8+0.3 


-10 


.68 


± 


0. 


.16 


4. 


.77+0. 


.88 


53 


.08+0. 


,25 


39 


.59 


-23 


,64 


-10 


.69 


6 


.2+2.9 


1308+277 


60+9 


-12 


.60 


± 


0. 


.18 


0. 


.52+0. 


.11 


52 


.88+0 


,31 


38 


.37 


-22 


,59 


-12 


.63 


12 


.6+7.2 


1319+428 


110+20 


-12 


.77 


± 


0. 


.16 


0. 


.43+0. 


.08 


52 


.03+0 


,30 


37 


.15 


-22. 


,79 


-12 


.90 


9 


.1+4.9 


1343+500 


4.0+0.5 


-10 


.71 


± 


0. 


16 


4. 


.61+0. 


,85 


52 


.76+0. 


,22 


39 


.17 


-23. 


,28 


-10 


.67 


7. 


.0+2.7 


1347+285 


21+4 


-12 


.31 


± 


0. 


13 


0. 


.73+0. 


,11 


51 


.13+0. 


,22 


36 


.35 


-23 


,85 


-12 


.72 


2 


.4+0.8 


1404+344 


2.8+0.2 


-10 


.59 


± 


0. 


16 


5, 


.23+0, 


,97 


53 


.00+0. 


.34 


39 


.72 


-23 


.15 


-10 


.50 


9 


.8+6.7 


1420+198 


43+6 


-12 


.27 


± 


0. 


.14 


0, 


.76+0. 


,12 


53 


.31+0. 


.39 


38 


.49 


-23 


.22 


-12 


.47 


4 


.4+3.3 


1441+262 


77+13 


-13 


.30 


± 


0. 


19 


0. 


.23+0. 


.05 


51 


.46+0 


.35 


36 


.49 


-23 


.40 


-13 


.46 


5 


.6+3.5 


1522+546 


43+7 


-12 


.33 


± 


0. 


.15 


0. 


.71+0. 


,12 


52 


.82+0. 


,30 


38 


.07 


-23 


,13 


-12 


.49 


5 


.1+2.7 


1533+557 


7.3+0.6 


-11 


.00 


± 


0. 


.16 


3 


,28+0, 


,60 


53 


.62+0. 


,30 


39 


.82 


-23 


.01 


-10 


.99 


7 


.9+4.7 


1547+215 


3.4+1.0 


-10 


.08 


± 


0. 


.18 


9. 


,44+1. 


.89 


53 


.02+0. 


,40 


39 


.45 


-22. 


,62 


-10 


.19 


5 


.8+4.1 


1549+628 


5.5+1.5 


-10 


.55 


± 


0. 


.15 


5. 


,50+0. 


.90 


52 


.91+0. 


,47 


39 


.34 


-22. 


,47 


-10 


.51 


9 


.3+4.8 


1609+660 


20+3.5 


-11 


.39 


± 


0. 


.18 


2. 


.10+0. 


.43 


53 


.37+0. 


.31 


39 


.35 


-22 


,19 


-11 


.36 


12 


.9+7.1 


1609+312 


12+3 


-12 


.08 


± 


0. 


.14 


0. 


.95+0. 


,15 


50 


.80+0, 


.23 


36 


.42 


-24. 


,00 


-12 


.39 


3 


.3+1.0 


1615+324 


47+9 


-12 


.24 


± 


0. 


15 


0. 


.79+0. 


,14 


52 


.25+0. 


.34 


37 


.78 


-22. 


,40 


-12 


.18 


10 


.8+6.3 


1618+177 


32+4 


-11 


.91 


± 


0. 


.15 


1. 


.15+0. 


,20 


53 


.49+0. 


.26 


38 


.83 


-23 


,10 


-12 


.07 


4 


.7+2.2 


1627+444 


24+6.5 


-11 


.35 


± 


0. 


13 


2. 


.20+0. 


,32 


53 


.23+0. 


.37 


39 


.03 


-22. 


.12 


-11 


.33 


9 


.6+4.0 


1658+302 


48+9 


-13 


.05 


± 


0. 


.17 


0. 


.31+0. 


,06 


50 


.89+0. 


.20 


36 


.10 


-23 


.67 


-13 


.12 


5 


.2+1.5 


1723+510 


20+3 


-11 


.54 


± 


0. 


18 


1, 


.76+0. 


,36 


53 


.48+0. 


.27 


39 


.55 


-22. 


.29 


-11 


.53 


16 


.0+7.7 


1726+318 


27+5 


-12 


.48 


± 


0. 


.14 


0. 


,60+0. 


,10 


52 


.75+0. 


.29 


38 


.14 


-23 


.76 


-12 


.67 


4 


.5+2.1 


1957+405 


13+3 


-10 


.39 


± 


0. 


.15 


6. 


,60+1. 


.18 


53 


.64+0 


.25 


39 


.39 


-22. 


.03 


-10 


.61 


4 


.7+1.6 


2019+098 


9.2+1.2 


-11 


.14 


± 


0. 


.10 


2. 


,80+0. 


,32 


53 


.33+0 


,27 


39 


.28 


-23 


,20 


-11 


.22 


5 


.6+2.7 


2104+763 


7.7+0.6 


-10 


.87 


± 


0. 


.14 


3 


,80+0. 


.63 


53 


.31+0 


,31 


39 


.06 


-23 


,28 


-11 


.25 


2 


.9+1.9 


2145+151 


6.4+1.1 


-10 


.75 


± 


0. 


13 


4 


,38+0. 


,64 


53 


.61+0. 


,38 


39 


.87 


-22 


,70 


-10 


.90 


7 


.4+3.9 
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Table 3: Observational and model 
cocoon 

Parameter 



parameters characterizing the source's 



Symbol Dimension 



Observational parameters from radio maps and spectra 

Apparent (projected) linear size D [kpc] 

axial ratio AR [dimensionless] 

observed volume V [kpc 3 ] 

1.4- GHz luminosity P 1A [WHz^sr- 1 ] 

source redshift z [dimensionless] 



Physical parameters derived directly from the above data 

equipartition magnetic field strength B eq [nT] 

equipartition energy density u eq [Jm~ 3 ] 

equipartition emitted energy U cq = u eq V [J] 



Physical parameters assumed in the model 

Central core radius ao [kpc] 

exponent of density profile j3 [dimensionless] 

exponent of particle energy distribution p [dimensionless] 

ratio of jet-head to cocoon pressure 7\c = Ph/Pc [dimensionless] 
adiabatic indices of ambient medium, 

cocoon, and magnetic field, respectively T x , r c , T B [dimensionless] 

source axis inclination 9 [°] 



Physical parameters fitted with the model for a given age 

Jet (constant) power Q [W] 

central core density po [kgm 



-31 



cocoon pressure p c [N m 

energy density u c [Jm" 

total emitted energy E tot = u c V c [J] 



- 2 ] 

-31 
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Table 4: The correlations between (log) D and t or Q or p when other 
parameters are kept constant 
Correlation r X y r XY /u Pxy/u 

T XY/V PxY/V 









TxY/W 


PxY/W 


T XV/TIVW 


PxV/TIVW 


D 


-t/Q 


+0.804 


+0.941 


<0.001 






D 


-t/po 




+0.805 


<0.001 






D 


-t/l+z 




+0.839 


<0.001 






D 


- t/Qo,p ,l+z 








+0.952 


<0.001 


D 


- Qo/Po 


-0.065 


+0.015 


0.898 






D 


- Qo/t 




+0.810 


<0.001 






D 


- Qo/l+z 




+0.450 


0.001 






D 


- Q /po,t,l+z 








+0.866 


<0.001 


D 


- Po/Qo 


-0.211 


-0.256 


0.031 






D 


- Po/t 




-0.080 


0.505 






D 


- Po/l+z 




-0.125 


0.298 






D 


- po/QoJA+z 








-0.716 


<0.001 



Table 5: The correlations between (log) P 1A and t or Q Q or 1+z when other 
parameters are kept constant 



Correlation 




txy/u 
rxv/v 


PxY/U 
PxY/V 


r XY/UV 


PxY/UV 


Pia — t/ Qo 
P 1A - t/l+z 
PiA-t/Q ,l+z 


-0.718 


-0.647 
-0.351 


0.001 
0.002 


-0.623 


0.005 


PlA — Qo/t 
PlA ~ Qo/l+Z 
PlA~Qo/t,l+Z 


+0.940 


+0.926 
+0.854 


<0.001 
<0.001 


+0.900 


<0.001 


Pia-(1+z)/Qo 

P 1A -(l+z)/t 

PiA-(l+z)/Q ,t 


+0.771 


+0.270 
+0.515 


0.023 
<0.001 


-0.153 


0.21 
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Table 6: The correlations between (log) AR and t or Q or p when the other 

parameters are kept constant 

Correlation r XY r X Y/u Pxy/u 

TXY/V PxY/V 









T XY /W 


PxY/W 


T XY/TJVW PxV/TJVW 


AR- 


t/Qo 


+0.328 


+0.615 


<0.001 




AR- 


t/po 




+0.445 


<0.001 




AR- 


t/l+z 




+0.488 


<0.001 




AR- 


t/Q ,po,l+z 








+0.513 <0.001 


AR- 


Qo/t 


+0.235 


+0.570 


<0.001 




AR- 


Qo/Po 




+0.121 


0.314 




AR- 


Qo/l+z 




+0.445 


<0.001 




AR- 


Qo/t,p ,l+z 








+0.421 <0.001 


AR- 


Po/t 


+0.225 


+0.380 


0.001 




AR- 


Po/Qo 




+0.182 


0.129 




AR- 


Po/l+z 




+0.301 


0.011 




AR- 


Po/t,Q ,l+z 








+0.209 0.085 



Table 7: The correlation between (log) p c and D or t or 1+z where the other 

parameters are kept constant 

Correlation r XY r X Y/u Pxy/u 

r XY/V PxY/V 
r XY/UV PxV/UV 

p c -D/t -0.792 -0.234 0.051 

p c -D/l+z -0.806 <0.001 

p c -D/t,l+z -0.533 <0.001 

p c -t/D -0.917 -0.772 <0.001 

p c - t/l+z -0.826 <0.001 

p c -t/D,l+z -0.467 0.001 

p c -(l+z)/D +0.785 +0.733 <0.001 

p c -{l+z)/t +0.145 0.23 

p c -(l+z)/D,t +0.324 0.005 
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